Quantum motion and level quantization of a skyrmion in a pinning potential in chiral magnets 
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A new topological excitation called skyrmion has been observed experimentally in chiral magnets without 
spatial inversion symmetry. The dynamics of a skyrmion is equivalent to an electron moving in a strong magnetic 
field. As a skyrmion involves large number of spins, it is not clear whether there exist observable quantum 
effects. In this work, we study the quantum dynamics of a skyrmion in a pinning potential. Without a pinning 
potential, the skyrmion performs cyclotron motion due to the strong emergent magnetic field originating from 
the Berry phase of spins, and all skyrmions occupy the lowest Landau level. Their wave functions are strongly 
localized in a region with radius less than 1 A when no driving force is applied. Thus in most circumstances, 
the quantum effects of a skyrmion are weak. In the presence of a pinning potential, the lowest Landau level for 
skyrmions is split into quantized levels, characterized by the orbital moments. The transition between quantized 
levels can be observed experimentally by microwave absorption measurements in low temperature region. The 
quantum effects are more prominent for a skyrmion with a small size, which can be realized in magnets with a 
large Dzyaloshinskii-Moriya interaction. 

PACS numbers: 75.25.-j, 75.70.Kw, 75.78.-n, 75.85. +t 



Recently stable triangular lattice of skyrmions in chi- 
ral magnets without inversion symmetry has been observed 
in conductors)^ such as MnSi, Feo.5Coo.5Si and also in 
insulator, 7 Cu20Se03, by the neutron scattering, Lorentz 
transmission electron microscope and spin-polarized scan- 
ning tunneling microscopy. In the temperature-magnetic field 
phase diagram, there are magnetic spiral phase, skyrmion lat- 
tice, ferromagnetic and paramagnetic states. In bulk sam- 
ples, the skyrmion phase occupies a small portion of the phase 
diagram, 1 while in thin films skyrmion lattice occupies a large 
area in the phase diagram. 5 For their promising applications 
in spintronics, skyrmions have attracted considerable atten- 
tion since the experimental observations. 

A skyrmion is a mesoscopic topological excitation in the 
ferromagnetic ground state. In skyrmion the spins wrap a 
sphere when one goes from the center of the skyrmion to 
the infinity, see Fig. [T] (a). When an external magnetic 
field is applied, skyrmions can be stabilized either by the 
Dzyaloshinskii-Moriya (DM) interaction 8 10 or by the dipole- 
dipole interaction 11 ! 12 !. Due to the topological nature, the 
skyrmions inside samples are stable against weak perturba- 
tions and disorder. Typically the size of skyrmion is about 10 
to 100 nm and involves 10 3 to 10 5 spins. In the conducting 
compounds where the Hund coupling is much stronger than 
the Fermi energy, the skyrmion can be driven by a spin polar- 
ized electric current through the spin transfer torque^E] i e 
transferring of the conduction electron's magnetic moments 
to localized spin moments; For the insulator Cu20Se03 with 
a strong magnetoelectric coupling, 16 it was proposed to drive 
the skyrmion by an electric field gradient 17 because skyrmion 
carries electric polarization. It was shown experimentally that 
the pining of skyrmion by quenched disorders is extremely 
weak. 18 20 The depinning current density for skyrmions is 
found to be about 10 5 ~ 10 6 A/m 2 , four order of magnitude 
smaller than that for the magnetic domain walls. 

Since a skyrmion involves huge amount of spin degrees of 
freedom, at the first sight one might expect that the quan- 
tum effects are weak. To date, most theoretical work on 



skyrmion treats it as a classical topological excitation. If 
these spins change collectively, observable quantum behav- 
ior of skyrmion may arise under certain conditions, and it is 
not clear how to observe the quantum effects in experiments. 
In the present work, we will show that the quantum motion of 
skyrmion in the presence of defects can result in energy level 
quantization. The quantized levels can be observed experi- 
mentally by microwave absorption or by an ac driving force. 

For a skyrmion, besides the extended modes (continuum 
of the spin wave spectrum), additional localized modes inside 
the skyrmion (internal modes) have been identified. 21 One is 
the Goldstone mode corresponding to the translation motion 
of skyrmion. The other internal modes, such as the breathing 
mode, are gapped. If the internal modes are not excited, we 
can consider only the translational motion of skyrmion and 
treat skyrmion as a rigid particle. We consider the two dimen- 
sional motion of skyrmions in a thin film. In this case, the 
equation of motion for skyrmions is governed by 22 
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where v,- = (v, v , v, y ) is the skyrmion velocity. The term on 
the left-hand side is the dissipative force with a damping co- 
efficient 77. It accounts for the damping of skyrmion motion 
and is induced by the underlying damping of the spin relax- 
ation. The mass is zero for a rigid skyrmion. When the in- 
ternal modes are excited, for instance when a skyrmion passes 
through a defect 23 , there will be additional contributions to the 
damping and inertial term 2425 . In this case, we need to add a 
term m s d,\j to the left-hand side of Eq. ([1} with m s being the 
skyrmion mass. The first term at the right-hand side of Eq. ^ 
Fm = 4ny zX V/ is the Magnus force per unit length, which is 
perpendicular to the velocity. Here £ is a unit vector normal to 
the film, and y - o 3 /(M) with a the lattice constant of mag- 
netic unit cell of chiral magnets and S the spin of the localized 
moments. Further F ss is the pairwise interaction between two 
skyrmions at r, and Tj. The term F p is the interaction between 
a skyrmion at r,- and a quenched disorder at r,. The last term 



Frf is the driving force. 

In conductors the driving force is due to the spin-transfer 
torque exerted by a spin polarized electric current J. It is given 
by Frf = z x JOo/c which is nothing but the Lorentz force, 
arising from the emergent quantized magnetic flux <t>o = hc/e 
carried by the skyrmion in the presence of a finite current. For 
insulators, the force can be produced by the gradient of elec- 
tric field E and it is given by F^ - P • d p E, where P is the 
total electric polarization inside a skyrmion and /u — x, y!^ 
Equation ([T| also describes the dynamics of vortices in type 
II superconductors in the presence of magnetic fieldsP§"2UFor 
vortices, the dissipative force is dominant over the Magnus 
force in most parameter space; while for skyrmions, the dis- 
sipation is due to the damping of precession of spins and is 
usually weak, thus the Magnus force is dominant. 

In real films the defects in the crystal structure cause pin- 
ning of skyrmions. Pinning potential can also be introduced 
intentionally with its strength and size being controllable in 
experiments (for example creation of columnar defects). For 
a properly engineered pinning potential, the cyclotron motion 
of skyrmion in the pinning potential in the quantum region 
produces experimentally observable consequences. 

Equation ([T} also describes an electron moving in an ef- 
fective magnetic field H* - 2<E>oc//a 3 for S = 1. For film 
thickness of d = 10 nm, the field is H* = 10 6 T for a * 3 A. 
With such a huge emergent internal field H", skyrmions are 
strongly confined in their cyclotron orbital with radius £ - 
a -\Ja/{4nd) ~ 0.15 A. Meanwhile skyrmion density can be 
controlled by applied magnetic fields. Near the phase bound- 
ary between the ferromagnetic phase and the skyrmion lattice 
phase, the density of skyrmion is small. In these cases, we can 
neglect the interaction between skyrmions putting F ss - 0. 
Neglecting the weak dissipation and using the canonical quan- 
tization, the Hamiltonian corresponding to Eq. ([T} with an 
inertia term, m s d,\i, is 

Wr) = l - [(-id x + A x f + (-id, + A,) 2 ] + U p (r) - F d x, (2) 

where r = (x, y) and U p is the pinning potential. The driving 
force is assumed along the x direction. Here dimensionless 
units are used. Energy is in units of 1i(L) c = 47rh 2 d/(a 3 m s ); 
length is in units of (; magnetic field is in units of H* and 
force is in units of fuo c jl. Since the mass of skyrmion is due 
to the excitation of internal modes, we expect the mass of the 
skyrmion is small. In our calculations we will take m s — > 
limit at the final step. Without driving force, Fj - 0, and at 
pinning potential U p = 0, the energy levels for the Hamilto- 
nian in Eq. |2]i are the well-known Landau levels with spacing 
%cj c . In the limit m s — > 0, it is sufficient to consider the lowest 
Landau level (LLL), where the skyrmion wave function under 
the symmetric gauge A = (y, -x)/2 is 
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¥„(0,r) = 



: exp(m6» - r 2 /2), 
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with an integer n corresponding to the orbital momentum 
along the z axis. As shown in Fig. [T] (b), the skyrmion 
wave function is strongly localized in the region with a ra- 
dius £, which is just the zero-point motion of skyrmion. As £ 
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FIG. 1. (color online) (a) Schematic view of the spin texture in a 
skyrmion. (b) The dynamics of skyrmions is equivalent to electrons 
moving in a strong magnetic field. The skyrmion performs cyclotron 
motion and its wave function is strongly localized in the absence 
of an external drive, (c) In the presence of a pinning potential, the 
LLL of skyrmion is split into discrete levels. Under a microwave 
irradiation, the skyrmion can be excited to a high energy level. 



is much smaller than the skyrmion size, the quantum effects 
of skyrmion in most cases are weak. The density of states 
in the LLL is given by N - H*^/^ - 2N S with J{ being 
the lateral area of the film and N s the total number of spins, 
which is far more bigger than the number of skyrmions. Thus 
it is impossible to fill the LLL by skyrmions, and the integer 
quantum Hall effect for skyrmions is absent. Nevertheless, it 
might be possible to have the fractional quantum Hall effect if 
the interaction between skyrmions is taken into account. 

When a driving force is applied, Fd > 0, skyrmion wave 
function is delocalized and is given by 
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where k y - 2mn/L y with L y being the length along the y di- 
rection and m an integer. These states are delocalized along 
the strip with the center at x — Fd - ky with width £. The 
corresponding energies are Ef, v = k y F d - F^/2. 

In the presence of the pinning potential U p and without a 
driving force, F d - 0, the degeneracy of the LLL with differ- 
ent orbital moment n is lifted, see Fig. |T](c). The energy for 
the orbital moment n is 

2 f°° 
E„ = - drU p ( r )r 2 " +l exp(-r 2 ). (5) 
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For convenience of calculations, we assume the pinning po- 
tential has the form 



U p (r) 



-V, 



r 2 + r 2 ' 



(6) 



where Uq is the pinning strength and r p is the size of the pin- 
ning potential. For columnar pinning potential with the size 
of order of the skyrmion size as shown below (see Fig. |2), 
Uq ~ J 2 x d/(Da) and r p /a ~ J ex /D, where / ex is the ex- 
change interaction and D is the DM interaction. 29 Note that 
the columnar pins can be created by heavy ions bombardment 
similar to that for superconductors. 30 31 The quantum level for 
a skyrmion in this pinning potential becomes equidistant and 
is given by 



AU = E n+l - E„ * 2U ({/r p ) 2 



(7) 



For the pinning potential Eq. |6|, we estimate AU ~ D/(2n). 
For MnSi, we have J ex w 3 meV and D « 0.3 meV 2 ^, thus 
AU ~ 0.05 meV. The quantized levels can be observed when 
the temperature is lower than T <k 0.5 K. For a stronger DM 
interaction AU becomes larger, which is helpful for experi- 
mental observations. 

Here we discuss the transition of skyrmion between differ- 
ent discrete levels in the pinning potential. The skyrmion 
can be excited to the higher level by an ac driving force 
F - F ac cos(cjt). When the frequency matches the level spac- 
ing ha> - AU, there will be resonant excitation of skyrmion 
to the high levels. The transition probability is determined by 
the matrix element F„,„ + i, 



F„,„+i = F. dc (n\x\n+ 1) 



F. dc ^£(2n + l)U 
2" +2 V«K«+ 1)! 



(8) 



For AU ~ 0.05 meV, the corresponding angular frequency is 
about co ~ 100 GHz. Experimentally the transition between 
different quantum levels of skyrmion can be excited by apply- 
ing microwave irradiation. For insulators like Cu20Se03, the 
microwave interacts strongly with the skyrmion. While for 
conductors like MnSi, the film thickness should be much less 
than the skin depth. For thickness of tens of nanometers, the 
microwave can fully penetrate through the film and provides 
strong coupling with skyrmions. 

Skyrmions need to be confined in the pinning potential 
during the measurements. In reality, skyrmions can escape 
the pinning potential via thermal and/or quantum fluctuations. 
The thermally-assisted escaping of the skyrmion can be min- 
imized by reducing the temperature well below the pinning 
strength, i.e. k B T <k Uq, with k B the Boltzmann constant. The 
quantum tunneling rate of skyrmion from one pinning center 
to its nearest pinni ng ce nter with a separation q is exp(— 5 g ) 
with S q ~ q 2 d/a 3 !P^21 Thus the quantum tunneling rate of 
skyrmions is negligible provided g » a. 

For a system with randomly distributed pining centers with 
varying pinning strength Uq and pinning size r p , AU may 
change from pinning centers to pinning centers. On average in 
the sample, the quantization levels become quasicontinuum, 
which poses a grand challenge for experimental observations. 




FIG. 2. (color online) Pinning potential U p versus the separation r„ 
between the center of the defect and the center of the skyrmion. The 
center of the defect is at r = and the size of the defect is f rf . The 
defect is modeled by J' cy {r) — [Vex _ Jd exp(-r/£/)] with J d = 0.8/ cx 
and £ d /a = J sx /D. The external field is H = 0.7£) 2 // cx . 



To overcome this difficulty, one may use a nanodisk of chiral 
magnets, where only one skyrmion can be stabilized in this 
confined geometry. 33 The pinning center can be introduced 
through heavy ion bombardment at the center of the disk. 

Finally we discuss the interaction between a skyrmion and 
a defect. We note a microscopic modeling of the defect has 
not been achieved yet. In the presence of defects in the crystal 
structure, both J ex and D becomes spatially dependent. More- 
over the magnetostriction will also present, which increases 
the effective size of defects. To model the interaction be- 
tween skyrmions and defects, one has to consider these ef- 
fects, which is beyond the scope of the present study. As 
a simple estimate sufficient for our purpose, we take a phe- 
nomenological approach to model the defect by a spatially de- 
pendent exchange interaction. The interaction between spins 
is given by the HamiltoniarP ' 10 ! 34 ! 35 ! 



W ' = " " Z ( 7 ex(r)n; n ; + D fa x nj)) - - J] H • n, (9) 
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where summation is over the nearest neighbors (i, /), D = 
Dx + Dy is the in-plane DM interaction and spins are rep- 
resented by a unit vector n. Without defects J' ex (r) = / ex , 
the skyrmion phase is thermodynamically stable in the region 
0.2D 2 /J ex < H < 0.8Z) 2 /7 ex PEU We model the defect by tak- 
ing J' ex (r) = [J ex - J d exp(-r/&)] with J d and £ being the 
strength and size of the defect respectively. Note that even for 
an atomic sized defect £/ « a, the pinning potential extends 
over the whole size of the skyrmion. To give an estimate, 
we assume that the distortion of the skyrmion induced by the 
defect is weak. We first obtain the solution of a skyrmion 
without the defect. We then substitute the obtained n, into 
Eq. [9] and calculate the self-energy of the skyrmion E s by 
fixing the separation r s between the center of the skyrmion 
and the defect. 22 The self-energy is just the pinning potential 



Up - E s . The results of U p as a function of r s are shown in 
Fig. [2] From Fig. [2] we can see the strength of the pinning is 
about J^ x d/{Da), which is a fraction of the self-energy of the 
skyrmion. The size of the pinning potential is r p ja as J ex /D, 
which roughly is just the radius of skyrmion. For a larger D, 
the radius of the skyrmion is smaller thus the quantum effects 
become more pronounced. 

To summarize, we have studied the quantum motion of a 
skyrmion in a pinning potential. The zero point motion of 
skyrmion is the cyclotron motion with radius less than 1 A. 
The skyrmions occupy the lowest Landau level with their 
wave functions strongly localized in the absence of a driv- 
ing force. Thus in most circumstances, the quantum effect of 
a skyrmion is weak. In the presence of a pinning potential, 
the degenerate lowest Landau level split into discrete levels. 
We have shown that for a properly engineered pinning po- 



tential, these discrete levels can be observed experimentally 
through microwave absorption measurements at low temper- 
atures. The quantum effects for skyrmions are more evident 
with a small skyrmion size, which can be achieved in magnets 
with a strong Dzyaloshinskii-Moriya interaction. 
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